A Measurement of the K° Charge Radius and a CP Violating Asymmetry Together 
with a Search for CP Violating El Direct Photon Emission in the Rare Decay 

Kl — > n + n~e + e~ 
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Using the complete KTeV data set of 5241 candidate Kl — > 7r + 7r~e + e _ decays (in- 
cluding an estimated background of 204±14 events), we have measured the coupling 
gen. = 0.163 ± 0.014 (stat) ± 0.023 (syst) of the CP conserving charge radius process and 
from it determined a K° charge radius of {r% ) = (-0.077 ± 0.007(stat) ± 0.011(syst))/m 2 . We 
have also determined a first experimental upper limit of 0.04 (90% CL) for the ratio r|^4 of 
the coupling for the El direct photon emission process relative to the coupling for Ml direct 
photon emission process. We also report the measurement of \qmi\ including its associated vector 
form factor \g M l\{l + ( M *_ m ^+% MkE ^ ) wh ere \g M i\ = 1.11 ± 0.12 (stat) ± 0.08 (syst) and 

ai/ffl2 = (— 0.744±0.027 (stat)±0.032 (syst))GeP 2 /c 2 . In addition, a measurement of the manifestly 
CP violating asymmetry of magnitude (13.6 ± 1.4 (stat) ± 1.5 (syst))% in the CP and T odd angle 
<j> between the decay planes of the e + e~ and it + it~ pairs in the Kl center of mass system is reported. 

PACS numbers: 13.20.Eb, 13.25.Es, 13.40. Gp, 13.40.Hq 



The emission of a virtual photon in the rare decay 
Kl — > 7r + 7r~e + e - proceeds via three main processes: 
bremsstrahlung, direct photon emission, and the charge 
radius process. The bremsstrahlung process takes place 
via the CP violating decay of a Kl — > 7r + 7r~ fol- 
lowed by emission of an electric dipole (El) photon by 
bremsstrahlung from one of the 7r's. The direct emis- 
sion process involves either the CP conserving or CP vi- 
olating direct emission at the primary decay vertex of a 
magnetic dipole (Ml) or a electric dipole (El) photon 
respectively. The CP conserving charge radius process is 
the transformation of a Kl — * Ks by emission of a vir- 
tual photon in a J=0 transition (forbidden in real photon 
emission) followed by the CP conserving decay of the Ks 
into 7r + 7r~. The charge radius coupling is related to the 
charge radius of the neutral kaon since the virtual photon 
acts as a probe of the K° in a way similar to the virtual 



photon in K° scattering from an atomic electron. The 
El and Ml direct emission and charge radius couplings 
are gsi, 9mi and gc'R- The matrix elements Q for the 
bremsstrahlung, Ml, El and charge radius processes are 
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where p+, p_, fc+, k, P are the 7r + , tt , positron, 
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electron, virtual photon, and Kl four momenta. The 
#0,1 are the 1=0,1 ir + ir~ strong interaction phase shifts. 
?7_l is the coupling of the Kl ir + ir~ decay. 

The KTeV E799-II experiment at Fermi National Ac- 
celerator Laboratory previously reported the first ob- 
servation of the rare four body decay mode Kl — > 
7r + 7r~e + e - based on 1% of the KTeV data. We have 
also made an initial measurement Q based on 36% of the 
KTeV Kl — > 7r + 7r~e + e~ data of a CP-violating asym- 
metry in the variable sin <p cos <p (where <j> is the CP- and 
T-odd angle between the e + e~ and 7r + 7r~ planes in the 
Kl cms). In addition, the measurement of the Ml direct 
photon emission coupling \qmi | including a vector form 
factor was reported in Ref. [3j. In this paper we report 
a measurement of the charge radius of the K° obtained 
from the coupling \gcR\ of the charge radius process of 
the Kl — > 7r + 7r~e + e~ decay. We also determined an up- 
per limit for the El direct photon emission in this decay. 
Finally, we present the measurements of the Ml direct 
emission process coupling and its form factor and the 
CP violating asymmetry in sin <f> cos <f>. 

The Kl — > 7r + 7r _ e + e~ data were accumulated during 
the 1997 and 1999 runs of the KTeV E799-II experiment. 
Differences in running conditions and spectrometer con- 
figuration can be found in Ref. 4]. The total KTeV E799- 
II Kl — * 7r + 7r~e + e~ signal of 5241 events, including an 
estimated background of 204 ± 14 events, obtained af- 
ter the analysis cuts described below, is shown in the 
7r + 7r~e + e - mass plot of Fig. ^ Note that data has been 
separately plotted for sin <j) cos 4> > and sin tj> cos <j> < 0. 
The CP violating asymmetry can be seen directly in the 
mass plot in the differing sizes of the two mass peaks. 

The KTeV four track trigger Q selected 3.9 x 10 8 
events from the 97 and 99 runs. Candidate Kl 
ir + n~ e + e~ events were extracted from these triggers by 
requiring events to have four tracks that passed track 
quality cuts and had a common vertex with a good vertex 
X 2 - To be designated as e ± , two of the tracks were re- 
quired to have opposite charges and 0.95 < E/p < 1.05, 
where E was the energy deposited by the track in the 
calorimeter, and p was the momentum obtained from 
magnetic deflection. To be consistent with a tt^ 1 pair, 
the other two tracks were required to have E/p < 0.90 
and opposite charges. To reduce backgrounds arising 
from other types of Kl decays in which decay prod- 
ucts have been missed, the candidate 7r + 7r _ e + e _ were 
required to have transverse momentum P 2 of the four 
tracks relative to the direction of the Kl be less than 
0.6 x 10~ 4 GeV 2 /c 2 . This cut was 94% efficient for re- 
taining Kl —> 7r + 7r~e + e~ . 

The major background to the Kl — > 7r + 7r~e + e~ mode 
was Kl — > ^tt^tt^ where -k° d was a Dalitz decay, 
7T° — > 7e + e~, in which the photon was not observed 
in the Csl calorimeter or the photon vetos. To reduce 
this background, all Kl — > 7r + 7r~e + e~ candidate events 
were interpreted as Kl — ► 7r + 7r _ 7r^, decays. Under this 



assumption, the longitudinal momentum squared (P^) v o 
of the assumed 7r° can be calculated in the frame in which 
the momentum of 7r + 7r~ is transverse to the Kl direc- 
tion. (PD^o was greater than zero for Kl — > Tr+ir^TTp 
decays except for cases where finite detector resolution 
resulted in a (P|) ff o < 0. In contrast, most of the 
Kl — > 7r + 7r~e + e~ decays had (Pl)^o < 0. The require- 
ment that all ir+ir-e+e- had (PgV < -0.025 GeV 2 /c 2 
reduced the Kl — > vt+tt - ^ background under the Kl 
peak to 177 events while retaining 94% of the signal. 
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FIG. 1: M 1T + n - e + e - invariant mass for events passing all 
Kl — » 7r + 7r~e + e~ physics cuts. The superimposed Kl mass 
peaks for sin cj> cos (j> > (white histogram) and < (gray his- 
togram) directly demonstrate the large CP violating asym- 
metry. There is no asymmetry in the backgrounds to the 
two peaks as demonstrated by their complete overlap of the 
distributions outside the kaon peak region. 

A second significant background was due to S° — > 
An — > p7r~e + e~7 where the photon was missed and the 
proton was misidentified as a ir + . There were 22 events 
of background after all cuts due to this decay. All other 
backgrounds were relatively minor. The largest was due 
to Kl — * 7r 7r°7r° with ir° — » e + e~e + e~. This mode con- 
tributed approximately four events to the background af- 
ter cuts. In addition, a potentially large background due 
to Kl —> 7r + 7r~7 decays in which the photon converted in 
the material of the spectrometer producing an e + e~ pair 
was eliminated by requiring M e + e - > 2.0 MeV/c 2 . The 
M e + e - cut retained 95% of the Kl — > 7r + 7r~e + e~ events 
with only one event contributing to the background. 

The final requirement of the Kl — * n + n~ e + e~ events 
was that 492 MeV/c 2 < M^ ee < 504 MeV/c 2 . The 
magnitude of the background under the Kl peak was 
determined by a fit to the simulated background mass 
distribution to the wings of the signal region. From this 
fit, a Kl — * TT + 7r~e + e~ signal of 5037 events above a 
background of 204 ± 14 events was obtained. 

We analyzed the Kl —> it + it~ e + e~ decays in a likeli- 
hood fit that used the matrix elements (eq. 1) of the 
model of Q with additional radiative corrections ap- 
plied to the final state particles using the PHOTOS pro- 
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gram |6| . We also found it necessary to include a depen- 
dence on the virtual photon energy in the Ml virtual pho- 
ton emission coupling in order to obtain agreement with 
the virtual photon energy spectrum £L* = E e + + E e - of 
the data (Fig.|2f). The Ml coupling \qmi\ wa s modified 
by a form factor 
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similar to that used in Ref |7j to describe Kl — > ir + Tr~j. 
Here M p is the mass of the p meson (770 MeV/c 2 ) and 
the photon energy has been replaced by E e + + E e - . 

The likelihood of a given event (see eq. 3 below), 
based on the matrix elements p,(xi,d) of the model of 
Ref. 0, is a function of the five independent variables 
xf. 4>, 8 e + (the angle between the e + and the tt + tt~ di- 
rection in the e + e~ cms), 9 n + (the angle between the 
7r + and the e + e~ direction in the 7r + 7r~ cms), M„+„-, 
and M e + e - . In addition, it depends on the values of 
the fit parameters a: a\jai and <7mi|, \g^\ > \9cp\ 
and nominal values for other model parameters such as 
77+_ = (2.286 ± 0.017) x 10~ 3 and = 43.51° ± 0.06°. 
The measured strong interaction phase shifts of the tt + it~ 
were taken from Ref. The likelihood was calculated 
using the selected Kl — > ir + n~e + e~ data sample of Nl> 
events and a large Monte Carlo event sample Nmc gener- 
ated with nominal values of the fit parameters do , passed 
through the spectrometer and reconstructed, and then 
reweighted with a new set of fit parameters d using the 
matrix elements p(xi,a). The likelihood fit to the five 
independent variables is shown in Fig. [21 along with the 
fit to Ery* . The charge radius process contributes to the 
higher mass M ee (as shown in the insert of Fig. 2c) while 
the Ml direct emission is determined by the shape of the 
M w + T - spectrum. 

The likelihood function used to perform the fit is 
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The best fit values were a\/a2 = (—0.744 ± 
0.027(stat)) GeV 2 /c 2 , \g M i\ = 1.11 ± 0.12(stat), \g C p\ = 
0.163 ±0.014(stat) and < 0.028 (upper limit due to 

statistical uncertainty only). The correlation (p = 0.924) 
between a\/a,2 and \§mi\ has been taken into account in 
determining their errors. 

The distribution of the quantity sin </> cos <fi (given by 
(n ee x hmr) • z(n ee ■ n T7r ), where the h's are the unit 
normals to the ee and tttt planes and z is the unit vector 
in the mr direction in the Kl cms) is shown in Fig. EK- 
The asymmetry of the sin <f> cos <j) distribution 
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yields (23.8 ±1.4 (stat))% before acceptance correc- 
tions. Using the fit of the model of Ref. 1] to the 
data to determine the acceptance, an asymmetry inte- 
grated over the entire Kl — > 7r + 7r _ e + e~ phase space of 
(13.6 ± 1.4 (stat))% was obtained, the largest such CP 
violating effect yet observed in kaon decay. The inter- 
ference between the Ml direct emission process and the 
bremsstrahlung process generates the asymmetry in the 
sin 6 cos 6 distribution. 




FIG. 2: Likelihood fit to the five independent variables a) 
sincj>cos4>, b) M 7r+7r - , c) M e+e - (the dotted curve in the insert 
in this figure shows the deficit of e + e~ pairs at high M ee if 
the charge radius process is ignored. The x 2 /dof for the fit 
of the model to the data increases from 0.85 to 1.6 if the 
charge radius process is left out.), d) 9^+, and e) 6 e + of the 
Kl — > iT + ir~e + e~ decay; f) i5 7 * defined as E e + + im- 
possible sources of false asymmetries were considered 
including those due to backgrounds and asymmetries in 
the detector. To check for detector asymmetries, a sam- 
ple of 15x 10 6 Kl — > 7r^7r~7T£) decays, which are expected 
to have no <fi asymmetry and which have similar topology 
to Kl — > 7r + 7r~e + e _ except for the presence of an extra 
photon in the Csl, was examined, and an asymmetry of 
((3.3 ± 2.6) x 10~ 2 )% was measured. 

Systematic errors on 121/02, |<?Afi|, \9cr\ and y^~r due 
to several sources are shown in Table [|] below. As shown 
in Table [I] the dominant systematic error is due to the 
variation of the fitted parameters resulting from varying 
the physics cuts used to select the Kl 
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and Monte Carlo events and repeating the fit procedure. 
Some analysis cut variations significantly increased the 
level of backgrounds to the K L — > 7r+7r~e + e~ mass peak. 
These cuts were separated from the other physics cuts 
and are labled as "background" in Table |U Finally, input 

parameters to the Monte Carlo such as rj^ ,<&-| , and 

the strong interaction tt + tt~ phases shifts <5o,i that were 
not included in the liklihood fit were varied by ±1ct of 
their published values to determine the uncertainty in 
the fit parameters due to their uncertainties. The total 
systematic errors in ai/a 2 , |gjwi|i \9cp\ an d were 
obtained by adding the systematic errors in quadrature. 

The systematic errors in the <f> asymmetry due to sev- 
eral sources are given in Table |n] below. The physics 
cut variations and background systematics of the <f> an- 
gle asymmetry have been determined as discussed above. 
The 77^ , <f>^ and <5o.i systematics are obtained as be- 
fore using the Ha uncertainties in these parameters. Ad- 
ditional uncertainties of the asymmetry due to the one a 
uncertainties of the fitted parameters are also included. 
All systematic errors of Table ITT1 are added in quadrature 
to obtain the total systematic error. 

In conclusion, the KTeV collaboration measured a 
charge radius coupling \gcp\ = 0.163 ± 0.014(stat) ± 
0.023(syst) which has been used to obtain, in a novel 
way a K° charge radius of {r 2 K a) = —S\gcB.\/MK = 
(-0.077±0.007(stat)±0.011(syst))(/m 2 ), consistent with 
the previous measurements of the K° charge radius [9lll(i 
\vi\ obtained in K° electron scattering and from a similar 
analysis of the Kl —> ir + Ti~e + e~ mode by NA48 0. We 
also set a first experimental upper limit on the presence 
of El direct photon emission in the Kl — > ir + ir~e + e~ 
mode of J|^4 < 0.04 (90%CL) including systematic er- 
rors. In addition, the Ml photon emission coupling was 
measured to be |<?mi| = 1-11 ± 0.12 (stat) ± 0.08 (syst) 
plus a vector form factor as given in equation (2) with 
ai/a 2 = (-0.744 ±0.027 (stat) ± 0.032 (syst)) GeV 2 /c 2 . 
Using a\ja 2 and |<?mi|j an average |<7mi| over the range 
of E y * was calculated to be 0.74±0.04. Finally, we made 
a measurement of a large CP-violating asymmetry in the 
distribution of T-odd angle (j) m Kl ~ * 7r + 7r~e + e~ de- 
cays of (13.6 ± 1.4 (stat) ± 1.5 (syst))% consistent with 
the theoretically expected asymmetry of Refs. 0,01. This 
result is consistent with our original measurement [3j and 
a later measurement by NA48 |l2j. 
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Uncertainty in 


Source 


ai/a2 




1 1 

\gcR\ 


9E 1 1 
9M1 


Monte Carlo Statistics 


0.002 


0.01 


0.001 


0.001 


Choice of initial MC parameters 


0.005 


0.02 


0.001 


0.001 


Skewing from input MC values 


0.000 


0.028 


0.002 


0.010 


Physics cut variations 


0.022 


0.041 


0.021 


0.018 


Background 


0.022 


0.05 


0.01 


0.008 


77^ Uncertainty 


0.0001 


0.01 


0.002 


0.0002 


$-1 Uncertainty 


0.0003 


0.002 


0.0002 


0.0005 


<5o,i Uncertainty 


0.001 


0.004 


0.001 


0.0003 


Total Systematic Error 


0.032 


0.08 


0.023 


0.023 



TABLE I: Syst. errors of a\/a2, \qmi\, \gcn\ and 



Source 


A Asymmetry (%) 


Physics cut variations 


0.71 


Background 


0.30 


7jn — Uncertainty 


0.163 


$-1 Uncertainty 


0.111 


5o.\ Uncertainty 


0.325 


l9Bll 
9M 1 1 


0.326 




IflMi \,ai/a 2 


0.335 


\gcn\ 


0.335 


Total Systematic Error 


1.46 



TABLE II: CP violating <f> asymmetry systematic errors 
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